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Abstract 

A holographic model of d + id superconductors based on the action proposed by Benini, 
Herzog, and Yarom a rXiv: 1006.0731] is studied. This model has a charged spin two field in an 
AdS black hole spacetime. Working in the probe limit, the normalizable solution of the spin 
two field in the bulk gives rise to a d + id superconducting order parameter at the boundary 
of the AdS. We calculate the fermion spectral function in this superconducting background and 
confirm the existence of fermi arcs for non-vanishing Majorana couplings. By changing the 
relative strength 7 of the d and id condensations, the position and the size of the fermi arcs 
are changed. When 7 = 1, the spectrum becomes isotropic and the spectral function is s-wave 
like. By changing the fermion mass, the fermi momentum is changed. We also calculate the 
conductivity for these holographic d + id superconductors where time reversal symmetry has 
been broken spontaneously. A non-vanishing Hall conductivity is obtained even without an 
external magnetic field. 



I. INTRODUCTION 



The mechanism of the high temperature superconductors (SC's) is an unsolved mystery 
in physics [Ij. High temperature SC's are layered compounds with copper-oxygen planes 
and are doped Mott insulators with strong electronic correlations. The pairing symmetry 
is unconventional and there is a strong experimental evidence showing that it is largely 
d-wave [2]. It is speculated that the pairing between electrons is mediated via strong 
anti-ferromagnetic spin fluctuations in the system. But the problem is difficult due to the 
strong-coupling nature of the system. Although significant progress has been made in the 
last many years [3l H], alternative approaches may be valuable to tackle the problem. 

An interesting alternative approach is the holographic correspondence between a grav- 
itational theory and a quantum field theory, which first emerged under the framework 
of AdS/CFT correspondence [5H7]. This method has provided us a useful and compli- 
mentary framework to describe strong interacting systems without the sign problem (see 
e.g. [8HT5]). In the original top-down approach within the AdS/CFT framework, both 
the gravity side and the field theory side were precisely know in string theory and gives 
us much deeper insight on this correspondence. But later in bottom-up approach we as- 
sume that the correspondence exists among the different pair of theories and try to make 
predictions from one side of the correspondence. 

Recently, a gravitational model of hairy black holes [HI [17] has been used to model 
s-wave high temperature SC's |18H2T]. In those class of models the Abelian symmetry of a 
complex scalar field is spontaneously broken (i.e. the Higgs mechanism) below some crit- 
ical temperature. The Meissner effect was soon observed by including a magnetic field in 
the background |22l [23] . A similar construction of holographic multiband superconductor 
has also been considered recently in plj . The multiband is coming from the condensation 
of a fundamental scalar field multiplet under non-abelian gauge group. The effect of the 
superconducting condensate on the holographic fermi surface has further been studied by 
calculating fermionic spectral functions [2SH2Z]- Interestingly, the properties of spectral 
function appeared to have similar behavior to that found in the angle resolved photo- 
emission experiment. Analogously, holographic p-wave superconductors have also been 
proposed by coupling a SU{2) Yang-Mills field to the black hole background, where a vec- 



tor hair develops in the superconducting phase [28H3T] . Properties of fermionic spectral 
function has also been studied in p-wave superconducting background [32]. Also, there 
are attempts to build holographic d-wave SC's by spontaneously breaking the Abelian 
symmetry of a charged spin two field [33H3S] and study different properties of the sys- 
tems [371439] . In addition to the bottom up approach mentioned above, there are various 
top-down constructions of those condensed matter like systems by considering D-brane 
configurations in the AdS black hole background in the string theory framework [lO] . 

In this work we will be interested in studying more detailed properties of holographic 
d-wave superconductors using the action of [Hil ES] • The most challenging problem to con- 
struct the holographic model of d-wave SC's is that the consistent action for the charged, 
massive spin two field in a general curve background is still not known. Although the 
action used in [351 EH] has some interesting features, such as having the right degrees of 
freedom and being ghost-free, it could be unstable or acausal for general gauge field con- 
figurations. It is argued that these problems may be cured by adding higher dimensional 
operators [SU] . 

In this work, we will study properties of the holographic d+id SC's based on the action 
of [351 [36] • This is motivated by the observations of spontaneous breaking of time reversal 
invariance in, for example, the YBCO high temperature superconductor [¥T1 - H5] . The 
time reversal breaking is thought to occur due to a complex combination of the d-wave 
condensates: d^2_y2 + id^y. One of the possible interesting consequences of this complex 
d-wave condensation is that the system has a Hall conductivity even in the absence of 
an external magnetic field [l6l HT]. By using the holographic framework, in this note 
we are interested in calculating the fermion spectral function, and the normal and Hall 
conductivity in the holographic d + id SC's. 

Our paper is organized as follow: In section|TT]we will first introduce our model following 
[HS] . We consider the effective action for the charged spin two field in a four dimensional 
AdS-Schwarzschild black hole. In order for the system to be ghost free, we will consider the 



probe limit. Subsequently in section [HT] we will obtain the normalizable background spin 
two field configuration which leads to a dual d + id superconducting phase of a strongly 
coupled system living on the boundary of AdS. Once we identify our required dual complex 
d-wave condensate background, we can study various linear transport properties of the 



system by considering linear order fermionic and gauge field perturbations. In section IV 



we will study the fermionic spectral function and discuss various properties depending 



upon various parameters such as the Majorana couphng, mass of the fermion etc. In 



section |VI] we will study the conductivity of the superconducting background. As we have 
mentioned earlier, because of spontaneous violation of time reversal symmetry due to the 
background condensation, we found a non-vanishing optical Hall conductivity. In the final 
section we will conclude with some remarks and discuss future directions. 



II. THE MODEL 

As we have mentioned in the introduction, the model that we are going to present is in 
the similar spirit of the holographic dual construction of an s-wave superconductor in an 
AdS black hole spacetime. It is well known that the phases of d-wave superconductivity 
can be described by a low energy effective theory of charged spin two tensor field in the 
framework of Landau-Ginzburg theory. At low temperature, the background solution of 
this charged spin two field yields the d-wave order parameter for the high temperature 
superconductivity. Analogously one particular construction of a holographic model of d- 
wave superconductor has been considered first in [33| l3l] by taking a charged symmetric 
traceless spin-2 field in an AdS-black hole background. Subsequently, in spite of having 
long-standing technical and also conceptual problems, a more refined form of the charge 
spin-2 tensor field effective action has been proposed and studied in the context of pure 
d-wave superconductor [36j. Our main goal in this report is to study the properties of 
the more general d + id superconducting background in AdS/CFT framework. One of the 
theoretical motivations to consider these kinds of systems is to discover some universal 
low energy properties of general strongly coupled field theories. Keeping this motivation 
in mind, we start by considering the action in the bulk containing the gravity part and 
the matter part as 

S=^ Jd'xV^(^(R+^^+cJj, (1) 

where R is the Ricci scalar, the term gives a negative cosmological constant and L 
is the AdS radius which will be set to unity in the units that we use. = SttGn is the 
gravitational coupling. The Lagrangian for the matter fields includes a spin two field and 
a spin half fermion field. Both of them are charged under U(l) gauge field and can be 
massive. We write 

= Cb + Cf . (2) 



The consistent construction of a charged spin two field theory in curved spacetime 
background is a long-standing problem. There have been lot of studies in constructing 
consistent interacting higher spin field theories [IH]. As we have mentioned, the authors 
in Ref. have proposed a unique Lagrangian for a charged spin-2 field in AdS space 
with the motivation in studying its dual field theoretic properties. In this report we will 
adopt their construction with the Lagrangian 

= -\D,^^,\^ + 2|Z}^y.'^n' + \DM'' - [D^v^^'D,^ + c.c] 

- m2(|^^j2 _ |^|2) ^ 2i?^,,,^*^V^' - R,uV*'^Vl (3) 

where (p^^^ is the spin two field with yj^j, = (yj^^, ip = (p^^, is the covariant derivative 
{Da^^^ = {da + iqAa) ffj^iy in fiat space), d{= 3)is the spatial dimension of the bulk and 
F^i, is the field strength tensor of the gauge field. As it has been discussed in detail in 
the above action describes dynamics of the correct number of degrees of freedom 
only when the background is Einstein's manifold. So, the energy density of the higher 
spin field configuration should be very small compared to the background energy density. 
This essentially means we have to maintain the probe limit of spin-2 and fermion fields 
in a specific gravitation background. However, there still exists serious issues in dealing 
with the above action. For generic gauge field background, even though it has correct 
number of dynamical degrees of freedom, the equation motion loses its hyperbolicity or 
causality. The general belief is that these violations can be ameliorated by considering 
higher derivative operators in our Lagrangian. This requirement makes it very complicated 
to construct a fully consistent Lagrangian for the higher spin field in a curved background. 
However as has been argued in ^B], we will adopt the effective field theory point of view, 
where all these effects could be very small in the limit of very low gauge field strength. 
This is essentially the limit in which we will be considering in our subsequent discussions. 

In addition to above Lagrangian for the spin-2 field, we also consider the spin 1/2 
fermion Lagrangian as 

Cf = {T^'D^ -m^)^ + r^V^^^L'^D'^^ - 7]^T^'D'' (<^^^^") . (4) 

The bulk Gamma matrices F'^ satisfies the Clifford algebra {F^,F'^} = 2g'^'^. The U{1) 
gauge symmetry demands that the charges of the fermion and spin two field are related 



by = q. It is know that 77 dependent Majorana coupling helps the fermion spectral 
function to develop a gap. It has been argued that although there are more terms of the 
same dimension as these 77 terms, only 77 dependent terms gives rise to an anisotropic gap. 
Thus, we have dropped the other terms for simplicity. 
Under the field redefinition 

^' -)■ "if/q, 7] qr], (5) 

C^n can then be written as 

Crn Cm/q^- (6) 

Thus, as g — )■ oo while the Cm remains finite, we can work in the so called probe limit 
to ignore the back reaction of the the matter fields A^^, (/?^^, and ^. Also, in the probe 
limit, the q dependence of observables can be recovered through the above simple scaling. 
Thus, we will work in the probe limit from now on. In this limit, the gravitational field 
equation satisfies the Einstein equation 

2A 

where A = — This yields the AdS Schwarzschild black hole solution with the metric 

ds' = ^ [-f(z)dt^ + dxU + ^ j , (8) 

with 

where z = Zh is the black hole horizon and 2; = is the boundary. The Hawking temper- 
ature can be expressed as 

T^j^. (10) 

An important point to note is that, the action is invariant under the scaling transforma- 
tions 

{t,x,y,z,T) {ct,cx,cy,cz,T/c), 
{A^, ^) ^ (A^/c, cpjc\ cr^ , (11) 

which determine the conformal dimension of each field. In the following section, we will 
consider this black hole background in order to solve the equations of motion for the 
spin-2 and electromagnetic gauge field in the probe limit. 



III. d + id CONDENSATE 



As we know from the standard AdS/CFT dictionary, the bulk massive spin-2 field 
corresponds to a spin-2 operator in the boundary field theory with fixed conformal di- 
mension. The conformal dimension is fixed by the spin-2 mass. This boundary spin-2 
dual operator under boundary Lorentz transformations is to be identified with the d- 
wave order parameter. Our main interest here is to describe the d-wave SC in the dual 
boundary field theory. Now since the spacetime background that we have considered has 
a translation symmetry in the boundary direction, the condensation on the x-y plane on 
the boundary should also have translational invariance. However rotational symmetry 
should be spontaneously broken down to Z{2) with the d-wave like condensate changing 
its sign under a 7r/2 rotation on the x-y plane. To incorporate these features, we use an 
ansatz for the symmetric traceless 9?^^: 

^xx = -^yy = ^^iW' 

fxy = fyx = ^V'2(^), (12) 

and ip^j^ = ioY 11,1/ ^ x,y, and for the gauge field A^: 

A = A^dx^" = (p{z)dt. (13) 
Under a 6 angle rotation in the x-y plan. 
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(14) 



If ^xy/Vxx is ^6^1) then we can always rotate the coordinates in the x-y plan such that 
ip^y — and then make ip^^ real by a gauge transformation. Otherwise ip^^ and ip^y will 
both exist. The equations of motion of the gauge field and the spin two field are 

0" + —0' - 4l d^il' + </> = 0, (15a) 

Z Z J 

< + - —W. + - = 0, i = 1, 2. (15b) 

J z z^f 

The coupling to the fermion field through is dropped because there is no fermion 

condensate. 



Near the horizon {5z = z — Zh —> 0) , f{z) ~ —d5z/zh, and 

SzhSz 



which yields 

Cl 



2*^^' 



'4'i - C2,i + Cs^ilnSz. (17) 



The finiteness of ipj^ at the horizon sets c^^i = 0. 
Near the boundary z ^ 0, / — )■ 0, and 



Z2 



2 

< - - ^ 0. (18a) 



9 r 2 

z z^ 



has the asymptotic solution with two terms: 



One of the terms can be identified as the source and the other term identified as the 
corresponding condensate. In other words, the asymptotic behavior of the spin two field 
near the boundary is 

{0^) 



+ Z^ 



+ Oiz) 



(20) 



^2A - d 

where is the source, (Oi) is the condensate, and Oi is the field theory operator that 
tp^ couples to at the boundary, = A{A — d) and A is conformal dimension of (Oi) 
(because the conformal dimension for is zero). 

Depending on the size of m^, there could be different scenarios: 

(a) If > 0, A_ < 0, we need to set Cs^j = (the sourceless condition) to keep ipi 
finite at the boundary. Thus, (j) has the asymptotic solution 

0^/i + p2 + O(z2), (21) 



where the physical meaning of is the chemical potential while p is the corresponding 
charge density. We can always use the scaling of Eq. (11) to set /z = 1. In this case 



V'2(^)/V'i(^) is a z independent constant. Because the combination C4,2V'i ~C4,iV'2 vanishes 
near the boundary. It continues to vanish for all z by the equations of motion. 

(b) If > > —9/4, A-t > 0, one can either chose the or term to be the 
condensate. If we choose the term to be the condensate for iIj^^ and choose the z^- 



term to be the condensate for '?/'2, then Eq.(21) is still true. But then '?/'2(-2)/V'i(-2^) is no 
longer a constant. 



(c) If m} = 0, A_ = 0, if we choose the z^+ term to be the source, then Eq.(21) is no 
longer true. Hence, the physical meaning in this case is unclear. The other choice of the 
source goes back to case (a). 

In case (b), the two condensates have different dimensions at the boundary. This might 
imply competition between the two order parameters in the system. However, the scaling 
dimension of any primary operator of a conformal field theory has a unitarity bound 
|49] . As has been argued in [36j, this unitarity bound constrains mass of a spin-2 field as 

> 0. Thus, in this report, we will just focus on case (a), where 

\Qif + \Q2f = l,^{z)=r{z). (22) 



We have solved the equations of motion for the spin-2 field ifj and gauge field using 
standard shooting algorithm by demanding yU = 1 and the normalizability of ip near the 
asymptotic boundary. We found the critical temperature Tc = 4~- for a fixed chemical 
potential and charge, below which non-trivial bulk profile for the spin-2 field component 
exists. According to AdS/CFT correspondence, the d-wave condensation (O) in the dual 
boundary field theory is identified with the coefficient of the normalizable solution of ip. 
In Fig, 1 we show the condensate field {{O) = \J\{Oi)'^ + |(02)|^) disappears above the 
critical temperature T^. Once we know the background condensation, it would be in- 
teresting to see the fermionic response function which essentially captures the properties 
of the background condensation. On the other hand, the electromagnetic perturbation 
is essential to study the conductivity of this background. In what follows we will first 
study the fermionic spectral function and discussed about the gap structure of the un- 
derlying strongly coupled system depending upon the type of fermionic coupling with a 
background we choose. Using the scaling relations in Eqs.([5]) and (11), one can show that 
the combination {qOY^^jj,/ p is independent of q and /z. 
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FIG. 1: {qO)^/^iJL/ p vs. T/Tc is shown. The combination in the y-axis is independent of q and 
/i. 

IV. FERMION SPECTRAL FUNCTION 



The equation of motion for the bulk fermion field looks like 



where ^p^ = D^(p^^ and the covariant derivative on the spinor field is 



(23) 



(24) 



with oj the spin connection, V^'^ = rt'^F'"], and with the vielbein indices underlined. We 
have chosen the bulk Gamma matrices representation 



-ia2 

Z(J2 

cxi 
ai 



as 

^3 
-20-2 

ia2 



(25) 



It is convenient to rescale ( = {—g ■ g^^y/^^iJ which removes the spin connection com- 
pletely, and to work in the momentum space with the Fourier components: 



-iuit+ik-x /-{u},k) 



iuit—ik-x ^(—u,—k) 



[z). 



(26) 



-1 r\ 



We then decompose the four-component spinor to two two-component spinors: ( = 
{(1X2)^ such that Eq. (23) can be rewritten as 



/=( — UJ, — kx)* I • >.( — w, — fcx)* 








(27) 



where 



Dia) = Vg^cTsd, - + + gc^t) V^^f^2 + iKVg^cr,. (28) 

Without losing generahty, we set ky = while the condensate can have arbitrary orienta- 
tions. 

The spectral function is defined by the imaginary part of the trace of the retarded 
green's function 

A{u,k) = lm[TT{GR)]. (29) 

To compute A{oj, k), we follow the method developed in [501452] and applied in [26] which 
leads to the simple "flow equation" : 



where 



- {iG) \iK^A + i{ujB + q^AtC)^/^^ + 2Kr^{g^^flH] {iG) 



(30) 
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(31) 



\ 1 / 

and S = V'' . The initial condition which satisfies the ingoing wave boundary condition at 
the horizon of an AdS black hole would be 
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I 2^-0 • 



(32) 




FIG. 2: The appearance of the fermi arcs in the fermion spectral density A{uj = 0, kx/qfJ-, ky/qfi) 
for normal d-wave superconductor (if^y = 0) with m = 0, ttiq = 0, tj = 0.15q and T = 0.66Tc. 



oj/qjJ 0.00 r - 




(a)6i = 0° 









— 1 












i 


























































































































o 




















































































































































O 




















































1 






















































i 


1 









0.0 0.1 0.2 0.3 0.4 



{h)e = 30° 



0.04 












1 










































0.02 














































r 










0.00 










— 1— 






















































+ 




























1 






















1 
























-4 






"1 










0.02 












1 










>- 




















































L 






































O.M 














--j- 






















1 





















0.0 0.1 0.2 0.3 0.4 



{c)e = 45° 



FIG. 3: The spectral function A{uj/qfj,,k/qfi) shown for various angles 9 with respect to the 
X-axis for m? = 2,T = O.QQTc and rj = 0.15q. 

Using the scaling relations in Eqs.(|5]) and (11 ), one can show that the spectral function 
A{u/q^,k/qfi) is independent of q and fi, provided ri/q is also q independent. 



A. Fermi Arc 



As mentioned above, when ^^y/^xx is real, we can always rotate the coordinate in 
the x-y plan such that ip^y = and make ip^^ real by a gauge transformation. The 
corresponding fermion spectral function A{u = 0,k/qfi) for ip^y = 0, r] = 0.15g, m = 0, 
= 0, and T = 0.66Tc is shown in Fig. [2j 

This case was first studied in Ref . [35] • It was found that the spectral function at w = 
does not vanish in the nodal direction but vanishes (up to some small values at small 



k = |k|) in other directions. This suggests that the system is gapped except in the nodal 
direction for T < Ti. At higher T, the ungapped directions (the so called "fermi arc") 
increase and eventually the system becomes ungapped above T2 (with T2 < Tc). This 
is reflected in the growing angular size of fermi arc with non-zero spectral function at 
Ti < T < T2. In this region, a typical spectral function A{oj/q^, k/q^) for various angles 
6 with respect to the x-axis, is shown in Fig. [3j Up to the small structure at low k, one 
can see the gap-like structure at 0° and 30° and the vanishing of the gap at 45°. At 0°, 
the spectral width is significantly smaller than the gap, suggesting it is a fermi liquid. 
One can also identify the fermi momentum at k/qfi ^ 0.3. In this model, the angular 
dependence of the fermi momentum is not built in. It would be interesting to generalize 
this model to incorporate this feature in the future. 

At this point we want to note the appearance of another inner fermi arc in the fermionic 
spectra function for a higher value of spin two field mass. This is also clearly visible in 
Fig. |3j As we decrease the mass of the background spin two field, this inner fermi surface 
disappears but outer fermi arc still exist. It is interesting to study this in more details and 
possibly ascribe any physical significance to this. Apparently the existence of two fermi 
surfaces suggest that in the boundary field theory, we may have two different characteristic 
fermionic degrees of freedom. 

B. Dirac Mass Dependence 

While it is interesting to see the fermi arcs in this holographic model, however, experi- 
mentally, the fermi arcs in high temperature SC's happen at Ti = Tc < T < T2, i.e. fermi 
arcs happen in the pseudo gap phase but not in the SC phase 

With the goal of eliminating the fermi arc in the SC phase in mind, we study the dirac 
mass {m(^) effect to the spectral function. We found that while tends to increase the 
gap, it also makes the nodal points not gapless any more. 

In Figjlj it is shown that the fermi surface shrinks and height of the peak is reduced 
when we turn on m^. If is sufficiently large, the spectral function will be gapped in 
all directions. 
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FIG. 4: The dirac mass dependence of the spectral function A{uj = 0, k/g/i). The fermi momen- 
tum decreases when is reduced, = 0, T = 0.66Tc and r] = 0.05q. 

V. d + id SUPERCONDUCTORS 

As mentioned above, if ^Pxy/Vxx is not real, then we cannot rotate the coordinates 
such that ip^y = 0, such that both ip^^ and ip^y will exist. A particular interesting case 
is ^xyl^xx — =1=^, whose ratio is unchanged under the rotation in Eq. (14). This means 



the difference between y?^^ and ^p^y in all directions is just a common phase which can be 
gauged away without physical consequences. This implies the fermion spectral function 
will be the same in all directions. 

In Figjsj we plot Aiu = 0, k) for different 

7 = (33) 

^xx 

As expected, the 7 = case is just a 7r/4 rotation of the 7 — ?■ oo case. Both of them are 
purely d-wave SC's. In the d + id SC with 7 = 1, the spectral function looks like a s-wave 
SC. 



VI. CONDUCTIVITY 

In this section we will calculate the conductivity by turning on an electromagnetic 
perturbation as previously discussed. According to the standard AdS/CFT dictionary, 
the gauge field perturbation in the bulk will lead to a boundary current. The boundary 
value of the perturbed gauge field becomes a source for this boundary current. 

The conductivity tensor aij can be defined as 

Ji = aijEj, (34) 
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FIG. 5: A{uj = 0, k/g/x) for d + id superconductors. For 7 = oo and 7 = 0, the spectral functions 
reduce to the d-wave case. For, 7 = 1, the spectral function is isotropic and is s-wave like. 
m2 = 0, T = O.eGTc, T/ = 0.05g, = 

where Ei and Jj are the external electric field and the induced current, respectively, in the 
i-direction {i = x,y). In linear response theory, aij is a current-current correlator which 
can be schematically denoted as aij ~ {Q \ [Jj, Jj] \ Q), where the matrix element denotes an 
ensemble average. Under a tt/2 rotation along the z-axis (R), R~^JiR = €ijJj, where ejj 
is an anti-symmetric tensor, and assuming the ensemble average is governed by properties 
of the ground state which is in general a d+id condensation, we have R\Q) = — 



Then, aij ~ (fi \RR-^ [Ji, Jj] RR-^\n) = {n \R~^ [J„ Jj] R\ Vl) = {n\[eikJk, ejiJi]\n). This 



implies a^^^: = ayy and a^y = -(Jy^. 

Under a parity (P^.) transformation with respect to the x-axis, the condensates trans- 
form as 



If {O2) = 0, as the case which is always achievable in a d-wave SC, then \Q) = ± \Q). 



Thus, a,, ~ {n\P,P-'[J„J,]P,P-^Q) = {n\P-'[J„J,]P,\n) ~ (-1)^-+^-(t,,. TWs 



yields axy = (Tyx = 0. In the case of d + id, however, \Q) is not an eigenstate of Px- Then, 




(35) 



in general axy = —Cyx 7^ 0. A similar argument using time reversal symmetry yields the 
same conclusion. So unlike the normal holographic d-wave SC, in addition to standard 
a the d + id SC could have a non- vanishing Hall conductivity axy component. 

In order to make our analysis simple, using the above symmetry we can diagonalize 



our conductivity tensor such that Eq,34 becomes 



(36) 



where 



J4 



Jx ± iJy ; E± 



Ex ± iEy ; a± = axx ± 



xy 



(37) 



As we have mentioned before, in order to calculate the conductivity, we need to consider 
the linearized perturbations in the bulk black hole background [53]. A consistent set of 
gauge field and spin-2 field perturbations are 
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^ty\" I ^ry\ 

where u is the frequency of the perturbation. Similarly, we can define the perturbation 
for the complex conjugate field ^v'^,^ also. In the above perturbation ansatz we have 
only considered zero momentum modes. We only limit ourself to consider the case where 

7 = Wxyl^xx = Q2/Q1 is real. 

With this perturbation ansatz we would like to obtain the linear order perturbation 
equations. Apparently as one can observe from [36] that when the background has only d- 
wave condensation, the equation of motion for two independent gauge field perturbations 
ax and ay are not mixed each other. From this one can arrive at the conclusion that there 
is no Hall current in the dual superfiuid phase. But in a. d + id SC, ax and ay couple 
which leads to Hall conductivity. The equation of motion in this general background 
is complicated to solve. But two decoupled set of equations emerge if we perform the 



-1 r> 



following field redefinition: 



Ctrn ~t~ ICLy-; dp CLx ^'^y^ 

^mt ~ ^ty ~~ "^^txi ^pt ~ ^ty "I" ^^tx^ 
^mz ^zy ^^zxi ^pz ^zy ~^ ^^zxi 

Smt ~ Sty ~ Ktxy Spt — Sty + ^Ste' 

^mz ^2j/ ^^zx' ^pz ^zy ~^ ^^zx- ("^^) 

Under this field redefinition one set of decoupled equations of motion become 

«m + y«m + + 1^ \{Ql ~ Ql)Cmt " {Q2 + Ql)imt] ' ^ ^ [(^2 " QljCLz 

- {Q2 + gi)C] + ^Yf^^'f - m [iQ2 - Ql)Cmz - iQ2 + Ql)Uz] = 0, 

imt + -^mt imt + ^^('^a + Qi)[^ + 2g(/>)a^ 

+ ^[(2(u; + g0)Cz + #'Cz]=O, 

z 

- iz^{uj + qcfy^'rat - «-(4w + 4g0 + qz^')^^t = 0, 

Cmt + -imt T-r imt + TTy(v2 " <5l)(-W + 2g0)a„ 



i 

~ 2 



(2(-u; + q4>)C + 



mz 



+ iz\-uj + q(t))C^ + 2|(-4a; + 4g0 + gz0')Ct = 0, (40) 

As we can see from the above set of equations, two equations for and Cmz algebraic 
in nature so we can use the third and the fifth equations to substitute Cmz 
other equations. Therefore, we have only three coupled fields equations for am,^rnt 
Crat to be solved. Similarly we have another set of equations for a-p, ipt^ Cpt^ ipz Cpz 
components. These equations can be obtained simply by replacing m — )■ p and Qi — )■ —Qi 
in the above. 

Now we will solve the above set of equations numerically by integrating them from 
the horizon to the boundary of the bulk AdS black hole spacetime. There exists a well 
defined procedure to calculate the retarded Greens function of the dual boundary field 
theory [511 ES]. At the horizon, the ingoing wave boundary condition of all the fields are 
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FIG. 6: The real part of the normal conductivity Re[(Txx] vs. T/Tc for holographic d-wave 
superconductors (7 = 0) with = 4. There are 5{u:) type supercurrent contributions in these 
curves that cannot be seen clearly. The hall conductivity axy vanishes. 

imposed to ensure causality: 

am = /"'^ ao H 

^mt — f ^ ^mtO + ■ ■ ■ ; ^mt ~ f ^ ^mtO + ' ' ' 

f = f~*t~^ f 4- ■ ■ ■ ■ f* = f~*l~^ (AD 

Where Oq, Cmto .^^^q are some arbitrary complex constants yet to be determined. Since 
the system in Eq. |40] is linear and homogeneous, we can set ao = 1 because this only 
affects the overall normalization of C,mto ^mto- effect on conductivity. So by 

demanding C,mt ^^"^ ^mt be normalizable near the boundary as was done in Ref. j36j . 
^^(0 Cmto uniquely fixed. Therefore we can solve for the conductivity. This can 
all be done without the shooting method. Because the system is linear and homogeneous, 
the problem is reduced to taking the correct linear combination of three independent 
solutions to eliminate the non-normalizable terms of and C,mf Once this is done, we 
can use the standard AdS/CFT technique to extract the retarded Greens function from 
the solution of the electromagnetic perturbation am- 
Near the boundary, has the asymptotic behavior 

am{z) = + al^z H , (42) 

where corresponds to an applied electric field in the boundary theory and corre- 





FIG. 7: Re\ 

^ xx\ of holographic d -\- id superconductors [m? — 4) for ^ — 1 with various T (left 
panel) and for T/Tc = 0.44 with various 7 (right panel). There are S{uj) type supercurrent 
contributions in these curves that cannot be seen clearly. 
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FIG. 8: The real and imaginary parts of the Hall conductivity, Re[axy] and Im[axy], for d + id 
superconductors with m? = 4, 7 = 1 and various T. 

sponds to the induced current. The standard definition of the Greens function [5ll |53j 
is 

cr+ = (Txx + icTxy = lim ^. (43) 

m 

Then the standard conductivity tensor component would simply be 

1 —i 
(Txx = + cr-), Gxy = Y^^+ ~ 

In the d-wave case with ip^y = 7 = 0, axy = 0, and we have reproduced the axx result 
of Refs. [5^ ES] in Fig|6} For 7 7^ 0, we have only shown the < 7 < 1 result in FigsjTjjo] 
because 7 — )■ I/7 under a 7r/2 rotation and then 

(^ij (7) = (^ij (1/7) • (45) 
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FIG. 9: Re[axy] and Imlaxy] for various 7 at = 4 and T/Tc = 0.44. Tc is independent of 7. 
axy = when 7 = 0. 

In each case, the real part of axx, ^Qio'xx), has a delta function type contribution S{u) from 
the super current if the system is in the superconducting phase. Also, Re{axx) is always 
non-negative. Below some frequency Ug, Re{axx) could have several peaks, suggesting 
there are several spin one resonance states with masses set by the scale qfi. Above 0;^, 
R^icxx) reaches its asymptotic value which is one. Numerically, Ug ~ qfi- Unlike the 
s-wave SC's, the fermi surfaces are not gapped everywhere for the d-wave or d + id SC's. 
These conducting electrons can respond to an electricity field of any frequency. Thus, 
there should be no gap in the conductivity. Our result also has this feature. 

For Re{axy), there is no 6{co) contribution and it does not have a definite sign. At the 
peaks of Re{axx), Reio-xy) vanishes (see Figs. [7|and[8|. (But Re{axy) vanishes does not 
imply a peak at Re{axx)-) This is because suppose one of the diagonalized conductivity 



cr+ is much bigger than the other in size, then according to Eq. (44), Re{axy) is nothing 
but lm{axx)- So it relates to Re{axx) as described above. 

In Figs. |8] and |9| the axy for various T and 7 is shown. It is interesting that within 
a small frequency window, the Hall conductivity Re{axy) could be very sensitive to the 
external electric field frequency with the Hall conductivity changes dramatically from a 
large negative value to a large positive value. Thus, it could be a good electromagnetic 
wave frequency sensor. About the 7 dependence, there is a nearly universal Re{axx) in 
Fig. [7] when u/qfi ~ 0.18. Also, the shape of axy in Fig. |9]is insensitive to 7. Those are 
curious features of this type of superconductors. 



VII. CONCLUSION 



We studied a holographic model of d+id superconductors based on the action proposed 
by Benini, Herzog, and Yarom [35] . The model contains a charged spin two field in an 
AdS black hole spacetime. Working in the probe limit, the normalizable solution of the 
spin two field in the bulk gives rise to a. d + id superconducting order parameter at the 
boundary of the AdS. 

We have calculated the fermion spectral function in this superconducting background 
and confirmed the existence of fermi arcs for non-vanishing Majorana couplings. Depend- 
ing on the relative strength 7 of the d and id condensations, we found that the position 
and the size of the fermi arcs are changing. Specifically when we take 7 = 1, the spectral 
functions become isotropic and are s-wave like. We also studied fermion mass effect. By 
changing the fermion mass, we saw the fermi momentum is changing. We have also cal- 
culated the conductivity for these holographic d + id superconductors where time reversal 
symmetry has been broken spontaneously. A non-vanishing Hall conductivity has been 
obtained even without an external magnetic field. 

As we know in a real high temperature superconductor, the fermi arc has been observed 
in a pseudo gap phase before the superconducting phase transition. So far in the existing 
holographic models of d-wave superconductor the fermi arc appears in the superconducting 
phase. This is an interesting open problem to construct such a holographic model. 
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